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ABSTRACT 



The results of different analyses of the dynamical parameters of the Universe are con- 
verging towards agreement. Remaining disagreements reflect systematic errors coming 
either from the observations or from differences in the methods of analysis. Compiling 
the most precise parameter values with our estimates of such systematic errors added, 
we find the following best values: the baryonic density parameter Qb h'^ = 0.019 ± 0.02, 
the density parameter of the matter component Qm = 0.29 ± 0.06, the density parameter 
of the cosmological constant Qx = 0.71 ± 0.07, the spectral index of scalar fluctuations 
Us = 1.02 ± 0.08, the equation of state of the cosmological constant wx < —0.86, and the 
deceleration parameter go = —0.56 ± 0.04. We do not modify the published best values 
of the Hubble parameter Hq = 0.73 ± 0.07 and the total density parameter Qq 11qx)2- 



1 INTRODUCTION 

Our information on the dynamical parameters of the Universe describing the cosmic ex- 
pansion comes from three different epochs. The earliest is the Big Bang nucleosynthesis 
which occurred a little over 2 minutes after the Big Bang, and which left its imprint in 
the abundances of the light elements affecting the baryonic density parameter Qf,- The 
discovery of anisotropic temperature fluctuations in the cosmic microwave background 
radiation at large angular scales (CMBR) by COBE-DMR |jl[, followed by small scale 
anisotropics measured in the balloon flights BOOMERANG and MAXIMA 0, by 
the radio telescopes Cosmic Background Imager (CBI) 0, Very Small Array (VSA) 
and Degree Angular Scale Interferometer (DASI) testify about the conditions in the 
Universe at the time of last scattering, about 350 000 years after Big Bang. The anal- 
yses of the CMBR power spectrum give information about every dynamical parameter, 
in particular Qq and its components Qb, and ^a, and the spectral index n^. For an 
extensive review of CMBR detectors and results, see Bersanelli et al. 0. Very recently, 
also the expected fluctuations in the CMBR polarization anisotropics has been observed 
by DASI i. 

The third epoch is the time of matter structures: galaxy clusters, galaxies and stars. 
Our view is limited to the redshifts we can observe which correspond to times of a few Gyr 
after Big Bang. This determines the Hubble constant, successfully done by the Hubble 
Space Telescope (HST) 0, and the difference Qx — in the dramatic supernova la 
observations by the High-z Supernova Search Team |TD[ and the Supernova Cosmology 



Project [|ri[]. The large scale structure (LSS) and its power spectrum has been studied in 
the SSRS2 and CfA2 galaxy surveys |jl2|, in the Las Campanas Redshift Survey |T^, in 
the Abell-ACO cluster survey |T^, in the IRAS PSCz Survey ||15[ and in the 2dF Galaxy 



2 



Redshift Survey [0,|T^. Various sets of CMBR data, supernova data and LSS data 
have been analyzed jointly. We shall only refer to global analyses of the now most recent 
CMBR power spectra and large scale distributions of galaxies. 

The list of other types of observations is really very long. To mention some, there have 
been observations on the gas fraction in X-ray clusters on X-ray cluster evolution 
T9[] , on the cluster mass function and the Lya forest [2^], on gravitational lensing 



on the Sunyaev-Zel'dovich effect ||22|, on classical double radio sources |2^, on galaxy 
peculiar velocities , on the evolution of galaxies and star creation versus the evolution 
of galaxy luminosity densities . 

In this review we shall cover briefly recent observations and results for the dynamical 
parameters iJo) ^b, ^m, ^x, ^o, n^, and go- In Section 2 these parameters are defined 
in their theoretical context, in Section 3 we turn to the Hubble parameter, and in Section 
4 to the baryonic density. The other parameters are discussed in Sections 5 and 6, which 
are organized according to observational method: supernova in Section 5, CMBR and 
LSS in Section 6. Section 7 summarizes our results. 



2 THEORY 

The currently accepted paradigm describing our homogeneous and isotropic Universe is 
based on the Robertson- Walker metric 



ds^ = c^dt' - d/2 = cMt^ - R(tf — ^ + a^de^ + a\in^e 



da 

and Einstein's covariant formula for the law of gravitation, 



G -^T (2) 

In Eq. (1) ds is the line element in four-dimensional spacetime, t is the time, R{t) is 
the cosmic scale, a is the comoving distance as measured by an observer who follows the 
expansion, k is the curvature parameter, c is the velocity of light, and 6, are comoving 
angular coordinates. In Eq. (2) G^u is the Einstein tensor describing the curved geometry 
of spacetime, T^jy is the energy-momentum tensor, and G is Newton's constant. 

From these equations one derives Friedmann's equations which can be put into the 
form 

R^ + kc^ SttG, 

= ^{pm+Px), (3) 

2R R^ + kc^ 8nG, 

^ + ^^ = -^(p™+Pa). (4) 
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Here p are energy densities, the subscripts m and A refer to matter and cosmological 
constant (or dark energy), respectively; and p\ arc the corresponding pressures of 
matter and dark energy, respectively. Using the expression for the critical density today, 

P.= 8^«S, (5) 

where Hq is the Hubble parameter at the present time, one can define density parameters 
for each energy component by 



O = p/pc • (6) 



The total density parameter is 



= + fir + • (7) 

In what follows we shall ignore the very small radiation density parameter VL^.. The matter 
density parameter Vt^ can further be divided into a cold dark matter (CDM) component 
ficDM, a baryonic component and a neutrino component Vt^. 

The pressure of matter is certainly very small, otherwise one would observe the galaxies 
having random motion similar to that of molecules in a gas under pressure. Thus one can 
set = in Eq. (4) to a good approximation. If the expansion is adiabatic so that the 
pressure of dark energy can be written in the form 

P\ = wxp\c^ , (8) 

and if dark energy and matter do not transform into one another, conservation of dark 
energy can be written 

Px + ?,Hpx{l + wx)^Q . (9) 
One further parameter is the deceleration parameter qQ, defined by 

RR R ^^^^ 



i?2 Rm ■ 

Eliminating R between Eqs. (4) and (10) one can see that go is not an independent 
parameter. 

The curvature parameter k in Eqs. (1), (3) and (4) describes the geometry of space: 
a spatially open universe is defined by A; = —1, a closed universe by A; = +1 and a flat 
universe by A; = 0. The curvature parameter is not an observable, but it is proportional 
to fio — 1, so if Q,Q is observed to be 1, the Universe is spatially flat. 
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3 THE HUBBLE PARAMETER 



From the definition of tlie Hubble parameter H = R/R one sees that it has the dimension 
of inverse time. Thus a characteristic time scale for the expansion of the Universe is the 
Hubble time 

TH = = 9.78h-^ X lOV- (11) 
Here h is the commonly used dimensionless quantity 

h = i/o/(100 km s"^ Mpc"^) . (12) 

The Hubble parameter also determines the size scale of the observable Universe. In 
time th, radiation travelling with the speed of light has reached the Hubble radius 

rn = thc = 3000/i-^Mpc. (13) 
Or, to put it differently, according to Hubble's non-relativistic law, 

z = Ho- , (14) 
c 

objects at this distance would be expected to attain the speed of light which is an absolute 
limit in the theory of special relativity. However, in special relativity the redshift z is 
infinite for objects at distance th receding with the speed of light and thus unphysical. 
Therefore no information can reach us from farther away, all radiation is redshifted to 
infinite wavelengths and no particle emitted within the Universe can exceed this distance. 

Our present knowledge of Hq comes from the Hubble Space Telescope (HST) Key 
Project 0. The goal of this project was to determine Hq by a Cepheid calibration of 
a number of independent, secondary distance indicators, including Type la supernovae, 
the Tully-Fisher relation, the fundamental plane for elliptical galaxies, surface brightness 
fluctuations, and type II supernovae. Here we shall restrict the discussion to the best 
absolute determinations of Hq, which are those from supernovae of type la. 

Visible bright supernova explosions are very brief events (one month) and very rare, 
historical records show that in our Galaxy they have occurred only every 300 years. The 
most recent one occurred in 1987 (code name SN1987A), not exactly in our Galaxy but in 
the nearby Large Magellanic Cloud (LMC). Since it now has become possible to observe 
supernova in very distant galaxies, one does not have to wait 300 years for the next one. 

The physical reason for this type of explosion (type SNII supernova) is the accumulation 
of Fe-group elements at the core of a massive red giant star of size 8-200 Mq which already 
has burned its hydrogen, helium and other light elements. Another type of explosion (type 
SNIa supernova) occurs when a degenerate dwarf star of CNO composition enters a stage 
of rapid nuclear burning to Fe-group elements. 
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The SNIa is the brightest and most homogeneous class of supernovse with hydrogen- 
poor spectra, their peak brightness can serve as remarkably precise standard candles 
visible from very far. Additional information is provided by the colour, the spectrum, and 
an empirical correlation observed between the time scale of the sharply rising light curve 
and the peak luminosity, which is followed by a gradual decline. Although supernovae are 
difficult to find, they can be used to determine Hq out to great distances, 500 Mpc or 
z ^ 0.1, and the internal precision of the method is very high. At greater distances one 
can still find supernovae, but Bubble's linear law (14) is then no longer valid, the velocity 
starts to accelerate. 

Supernovae of type II are fainter, and show a wider variation in luminosity. Thus they 
are not standard candles, but the time evolution of their expanding atmospheres provides 
an indirect distance indicator, useful out to some 200 Mpc. 

Two further methods to determine Ho make use of correlations between different galaxy 
properties. Spiral galaxies rotate, and there the TuUy-Fisher relation correlates total 
luminosity with maximum rotation velocity. This is currently the most commonly applied 
distance indicator, useful for measuring extragalactic distances out to about 150 Mpc. 
Elliptical galaxies do not rotate, they are found to occupy a "fundamental plane" in which 
an effective radius is tightly correlated with the surface brightness inside that radius and 
with the central velocity dispersion of the stars. In principle this method could be apphed 
out to z ~ 1, but in practice stellar evolution effects and the non-linearity of Bubble's 
law limit the method to z < 0.1, or about 400 Mpc. 

The resolution of individual stars within galaxies clearly depends on the distance to 
the galaxy. This method, called surface brightness fluctuations (SBF), is an indicator of 
relative distances to elliptical galaxies and some types of spirals. The internal precision 
of the method is very high, but it can be applied only out to about 70 Mpc. 

Observations from the HST combining all this methods M and independent SNIa 



observations from observatories on the ground [26| agree on a value 



ifo = 73 ± 2 ± 7 km s"^ Mpc'\ (15) 

Note that the second error in Eq. (15) which is systematical, is much bigger than the 
statistical error. This illustrates that there are many unknown effects which complicate 
the determination of Hq, and which in the past have made all determinations contro- 
versial. To give just one example, if there is dust on the sight line to a supernova, its 
light would be reddened and one would conclude that the recession velocity is higher than 
it in reality is. There are other methods such as weak lensing which do not suffer from 
this systematic error, but they have not yet reached a precision superior to that in Eq. (15). 
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4 THE BARYONIC DENSITY 

The ratio of baryons to photons or the baryon abundance is defined as 

r/ = ^ ~ 2.75 X 10"*^ rib h'^ (16) 

where A*";, is the number density of baryons and A^^ = 4.11 x 10^ is the number 
density of photons. Thus the primordial abundances of baryonic matter in the standard 
Big Bang nucleosynthesis scenario (BBN) is proportional to VL^h"^. Its value is obtained 
in direct measurements of the abundances of the light elements ^He, ^He, or D, ^Li 
and indirectly from CMBR observations and galaxy cluster observations. 

If the observed abundances are indeed of cosmological origin, they must not signifi- 
cantly be affected by later stellar processes. The helium isotopes ^He and '^He cannot 
be destroyed easily but they are continuously produced in stellar interiors. Some recent 
helium is blown off from supernova progenitors, but that fraction can be corrected for 
by observing the total abundance in hydrogen clouds of different age, and extrapolating 
it to time zero. The remainder is then primordial helium emanating from BBN. On the 
other hand, the deuterium abundance can only decrease, it is easily burned to ^He in 
later stellar events. The case of ''Li is complicated because some fraction is due to later 
galactic cosmic ray spallation products. 

Among the light elements the ^He abundance is easiest to observe, but also least 
sensitive to Q^h"^, its dependence is logarithmic, so that only very precise measurements 
are relevant. The best "laboratories" for measuring the ^He abundance are a class of low- 
luminosity dwarf galaxies called Blue Compact Dwarf (BCD) galaxies, which undergo an 
intense burst of star formation in a very compact region. The BCDs are among the most 
metal-deficient gas-rich galaxies known. Since their gas has not been processed during 
many generations of stars, it should approximate well the pristine primordial gas. 

Over the years the observations have yielded many confiicting results, but the data 
are now progressing towards a common value [^, in particular by the work of Yu. I. 



Izotov and his group. The analysis in their most recent paper |]2S|, based on the two most 
metal-deficient BCDs known, gives the result 

f)fe(^He) h"^ = 0.017 ± 0.005 {2ct CL) , (17) 

where the error is statistical only. Usually one quotes the ratio Yp of mass in ^He to 
total mass in and ^He, which in this case is 0.2452 with a systematic error in the 
positive direction estimated to be 2-4%. Because of the logarithmic dependence, this 
error translated to Qb h"^ could be considerable, of the order of 100% . 

The ^He isotope can be seen in the Milky Way interstellar medium and its abundance 
is a strong constraint on Qbh^- The ^He abundance has been determined from 14 years 
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of data by Balser et al. [^. More interestingly, Bania et al. |^| combined Milky Way 
data with the helium abundance in stars to find 



QbCRe) h"^ = 0.020 



+0.007 
-0.003 



{1(1 CL) 



'1^ 



There are actually three different errors in their analysis, and their quadratic sum gives 
the total error. The first error is from the observed emission-line that includes the errors in 
the Gaussian fits to the observed line parameters. The second error is from the standard 
deviation of the observed continuum data and the third error is the percent uncertainty 



of all models that have been used in the analyses of reference |29 
For a constraint on Qi, h"^ from ^Li, Coc et al. 



update the previous work of several 
groups. More importantly, they include NACRE data in their compilation, and the 
uncertainties are analysed in detail. There is some lack of information about the neutron- 
induced reaction in the NACRE compilation, but the main source of uncertainty for the 
lighter neutron-induced reaction (e.g. ^H(n, 7)^H and ^He(n, p)^H) is the neutron lifetime 
(for the present value see the Review of Particle Physics [0). However, there is no new 
information about the heavier neutron-induced reaction (e.g. ^Li) or for ^He(d, p)^He, 
but in this compilation the Gaussian errors have been opted from the polynomial fit of 



NoUett & Buries [35|. We quote Coc et al. 132| for 



fi;,(^Li) /i^ = 0.015 ± 0.003 (1(7 CL) . 



(19) 



The strongest constraint on the baryonic density comes from the primordial deuterium 
abundance. Deuterium is observed as a Lyman-a feature in the absorption spectra of 
high-redshift quasars. A recent analysis ||36| gives 



flbCR) h'^ = 0.020 ± 0.001 (la CL) , (20) 

more precisely than any other determination. Some systematic uncertainties remain in 
the calculations arising from the reaction cross sections. 
Very recently Chiappini et al. 



37] have redefined the production and destruction of 



^He in low and intermediate mass stars. They also propose a new model for the time 
evolution of deuterium in the Galaxy. Taken together, they conclude that Qi,h'^ > 0.017, 
in good agreement with the values in Eqs. (18) and (20). 

Let us now turn to the information from the cosmic microwave background radiation 
and from large scale structures. There are many analyses of joint CMBR data, in partic- 
ular three large compilations. Percival et al. combine the data from COBE-DMR [|T| 



MAXIMA m, BOOMERANG 0, DASI g], VSA § and CBI § with the 2dFGRS LSS 



data [^]. Wang et al.^l| combine the same CMBR data (except VSA) with 20 earlier 



CMBR power spectra, take their LSS power spectra from the IRAS PSCz survey |jT^, and 
include constraints from Lyman a forest spectra |^ and from the Hubble parameter 
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quoted in Eq. (15). Sievers et al. also use the same CMBR data as Percival et al.[ 
(except VSA), combine them with earher LSS data, and use the HST Hubble parameter 
quoted in Eq. (15) and the supernova data referred to in Section 5 as supplementary 
constraints. All these analyses are maximum likelihood fits based on frequentist statistics, 
so the use of the Bayesian term "prior" for constraint is a misnomer. 

Assuming that the initial seed fluctuations were adiabatic, Gaussian, and well described 
by power law spectra, the values of a large number of parameters are obtained by fitting the 
observed power spectrum. Here we shall only discuss results on Qh h"^ which is essentially 
measured by the relative magnitudes of the first and second acoustic peaks in the CMBR 
power spectrum, returning to this subject in more detail in Section 6. 

The data used in the three compilations are overlapping but not identical, and the 
central values show a spread over ±0.0003. This we treat as a systematic error to the 
straight unweighted average of the central values. Two compilations |^ consider 
models with and without a tensor component. Since the fits are equally good in both 
cases we take their difference, ±0.0008, to constitute another systematic error. We shall 
use this averaging prescription also in Section 6 to obtain values of other parameters. All 
the analyses can then be summarized by the value 



f^b(CMBR) h'^ = 0.022 ± 0.002 ± 0.001 (Icr CL) 
where the statistical error corresponds to references [E^, pi||. 



(21) 



Method 



7] 



Error 



References 



^He abundance 
^He abundance 
''Li abundance 
^H abundance 
CMBR ± 2dFGRS 



4.7 til X 

+2.2 
-1.2 ^ 



10-10 

5.4 X 10-1° 
5.0 X 10-1° 



0.017 ±0.005 

090 +° °°'^ 

U.UZU _0.003 

0.015 ±0.003 



5.6 ± 0.5 X 10-1° 0.020 ±0.001 



2(7 stat. only 
la stat. only 
la stat. only 
la stat.±syst. 



m 
m 

m 



0.022 ± 0.002 ± 0.001 la stat.± syst. [Ulill 



Table 1: The baryonic density parameter 

In Table 1 we summarize the results from Eqs. (17-21). From this table one can 
conclude that all determinations are consistent with the most precise one from deuterium 
||36|| . A weighted mean using the quoted errors yields 0.0194 ± 0.0008 which is dominated 
by deuterium. However, all light element abundance determinations generally suffer from 
the potential for systematic errors. As to CMBR, the statistical errors quoted in all 
compilations have been obtained by marginalizing, so they are certainly unrealistically 
small. We take a conservative approach and add a systematic error of ±0.002 linearly to 
each of the five data values before averaging. The weighted mean is then 
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= 0.019 ±0.002 , (22) 

in excellent agreement with all the uncorrected input values in Table 1. 

One further source of Qf, information is galaxy clusters which are composed of baryonic 
and non-baryonic matter. The baryonic matter takes the forms of hot gas emitting X- 
rays, stellar mass observed in visual light, and perhaps invisible baryonic dark matter of 
unknown composition. Let us denote the respective fractions fgas, fgai, and fbdm- Then 

fgas ± fgal fbdm "^Q ' i.'^^) 

where T describes the possible local enhancement or diminution of baryon matter density 
in a cluster compared to the universal baryon density. This relation could in principle be 
used to determine fib when one knows flm (or vice versa), since fgas and fgai can be mea- 
sured, albeit with large scatter, while fbdm can be assumed negligible. Cluster formation 
simulations give information on T |^,|^5| to a precision of about 10%. However, the 



precision obtained for Qb h'^ by adding several 10% errors in quadrature does not make 
this method competitive. 



5 SUPERNOVA la CONSTRAINTS 

In Section 3 we already mentioned briefly the physics of supernovae. The SN la observa- 
tions by the High-z Supernova Search Team (HSST) |]lOl and the Supernova Cosmology 
Project (SCP) ||rT| are well enough known not to require a detailed presentation here. 
The importance of these observations lies in that they determine approximately the lin- 
ear combination Vtx — VLm which is orthogonal to = + ^^a, see Figure 1. 

HSST use two quite distinct methods of light-curve fitting to determine the distance 
moduli of their 16 SNe la studied. Their luminosity distances are used to place constraints 
on six cosmological parameters: /i, ^a, (Jo, and the dynamical age of the Universe, to- 
The MLCS method involves statistical methods at a more refined level than the empirical 
template model. The distance moduli are found from a analysis using an empirical 
model containing four free parameters. The MLCS method and the template method give 
moduli which differ by about la. Once the distance moduli are known, the parameters h, 
ilmj ^A are determined by a maximum likelihood fit, and finally the Hubble parameter is 
integrated out. (The results are really independent of h.) One may perhaps be somewhat 
concerned about the assumption that each modulus is normally distributed. We have no 
reason to doubt that, but if the iterative analysis has yielded systematically skewed 
pdf's, then the maximum likelihood fit will amplify the skewness. 
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High-z Supernovae Search Team 
(Riess et al, 1998) 



Supernova Cosmology Project 
(Perlmutter et al., 1999) 



m 



Figure 1: The best fit confidence regions in {fim — ^a) plane in the analyses of the Supernova 
Cosmology Project (blue curves) |11] and the High Redshift Supernova Search Team (red curves) 
[10|. The diagonal line corresponds to a flat cosmology. Above the flat line the Universe is closed 
and below it is open. 
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The authors state that "the dominant source of statistical uncertainty is the extinction 
measurement". The main doubt raised about the SN la observations is the risk that 
(part of) the reddening of the SNe la could be caused by intervening dust rather than 
by the cosmological expansion, as we already noted after Eq. (15). Among the possible 
systematic errors investigated is also that associated with extinction. No systematic error 
is found to be important here, but for such a small sample of SNe la one can expect that 
the selection bias might be the largest problem. 

The authors do not express any view about which method should be considered more 
reliable, thus noting that "we must consider the difference between the cosmological con- 
straints reached from the two fitting methods to be a systematic uncertainty" . We shall 
come back to this question later. Here we would like to point out that if one corrects for 
the unphysical region < using the method of Feldman & Cousins EBf , the best value 



and the confidence contours will be shifted slightly towards higher values of f2o. This shift 
will be more important for the MLCS method than for the template method, because the 
former extends deeper into the unphysical Qm region. 

Let us now turn to SCP, which studied 42 SNe la. The MLCS method described 
above is basically repeated, but modified in many details for which we refer the reader 
to the source |11|. The distance moduli are again found from a analysis using an 
empirical model containing four free parameters, but this model is slightly different from 
the HSST treatment. The parameters and Qx are then determined by a maximum 
likelihood fit to four parameters, of which the parameters Ai b (an absolute magnitude) 
and a (the slope of the width-luminosity relation) are just ancillary variables which are 
integrated out (h does not enter at all). The likelihood contours in (Qm — ^x) plane of both 
supernovae projects (SCP and HSST) are shown in Figure-1. The authors then correct 
the resulting likelihood contours for the unphysical region < using the method of 
Feldman & Cousins Since the number of SNe la is here so much larger than in 

HSST, the effects of selection and of possible systematic errors can be investigated more 
thoroughly. SCP quotes a total possible systematic uncertainty to flf^""* and of 0.05. 

If we compare the observations along the line defining a flat Universe, SCP finds fl\ — 
= 0.44 ± 0.085 ± 0.05, whereas HSST finds Qx - = 0.36 ± 0.10 for the MLCS 
method and fix — = 0.68 ± 0.09 for the template method. Treating this difference as 
a systematic error of size ±0.16 the combined SCP result is 0.52 ± 0.10 ± 0.16. SCP and 
HSST then agree within their statistical errors - how well they agree cannot be established 
since they are not completely independent. We choose to quote a combined HSST and 
SCP value 

nx-nm = o.5±o.i , (22) 

which excludes a flat de Sitter universe with Qx — = 1 by 5a, and excludes a fiat 
Einstein - de Sitter universe with Qx ~ = — 1 by lOa. 
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6 CMBR AND LSS CONSTRAINTS 

The most important source of information on the cosmological parameters are the anisotro- 
pies observed in the CMBR temperature and polarization maps over the sky. The tem- 
perature angular power spectrum has been measured and analyzed since 1992 [|l], whereas 
the polarization spectrum is very recent [§] and has not yet been analyzed to obtain val- 
ues for the dynamical parameters. Given the temperature angular power spectrum, the 
polarization spectrum is predicted with essentially no free parameters. At the moment 
one can say that the temperature angular power spectrum supports the current model 
of the Universe as defined by the dynamical parameters obtained from the temperature 
angular power spectrum. 

Temperature fluctuations in the CMBR around a mean temperature in a direction a 
on the sky can be analyzed in terms of the autocorrelation function C(^) which measures 
the average product of temperatures in two directions separated by an angle ^, 

CW = (^(«)^(« + ^)) ■ (23) 

For small angles (ff) the temperature autocorrelation function can be expressed as a sum 
of Legendre polynomials Pi{&) of order the wave number, with coefficients or powers 

1 oo 

C{e) = — ^a£2(2£ + l)P^(cos^) . (24) 

1=2 

All analyses start with the quadrupole mode 1 = 1 because the i = monopole mode 
is just the mean temperature over the observed part of the sky, and the i = 1 mode is 
the dipole anisotropy due to the motion of Earth relative to the CMBR. In the analysis 
the powers aj are adjusted to give a best fit of C{9) to the observed temperature. The 
resulting distribution of aj values versus i is the power spectrum of the fiuctuations, see 
Figure 2. The higher the angular resolution, the more terms of high i must be included. 

The exact form of the power spectrum is very dependent on assumptions about the 
matter content of the Universe. It can be parametrized by the vacuum density parameter 
Qk = 1 — ^0; the total density parameter Qq with its components flm, ^x, and the matter 
density parameter Q^n withits components Q^, ^cdm, ^u- Further parameters are the 
Hubble parameter h, the tilt of scalar fiuctuations Ug, the CMBR quadrupole normal- 
ization for scalar fiuctuations Q, the tilt of tensor fiuctuations rit, the CMB quadrupole 
normalization for tensor fiuctuations r, and the optical depth parameter r. Among these 
parameters, really only about six have an infiuence on the fit. 
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In Section 4 we already noted that the relative magnitudes of the first and second 
acoustic peaks are sensitive to ilf,- The position of the first acoustic peak in multipole i 
- space is sensitive to Qq, which makes the CMBR information complementary (and in 
Qm, f^A - space orthogonal) to the supernova information. A decrease in corresponds 
to a decrease in curvature and a shift of the power spectrum towards high multipoles. An 
increase in fix (in Sat space) and a decrease in h (keeping fl^ h'^ fixed) both boost the 
peaks and change their location mi- space. 

Let us now turn to the distribution of matter in the Universe which can, to some 
approximation, be described by the hydrodynamics of a viscous, non-static fluid. In 
such a medium there naturally appear random fluctuations around the mean density p{t) , 
manifested by compressions in some regions and rarefactions in other regions. An ordinary 
fluid is dominated by the material pressure, but in the fluid of our Universe three effects 
are competing: radiation pressure, gravitational attraction and density dilution due to 
the Hubble flow. This makes the physics different from ordinary hydrodynamics, regions 
of overdensity are gravitationally amplified and may, if time permits, grow into large 
inhomogenities, depleting adjacent regions of underdensity. 

Two complementary techniques are available for theoretical modelling of galaxy for- 
mation and evolution: numerical simulations and semi-analytic modelling. The strategy 
in both cases is to calculate how density perturbations emerging from the Big Bang turn 
into visible galaxies. This requires following through a number of processes: the growth 
of dark matter halos by accretion and mergers, the dynamics of cooling gas, the trans- 
formation of cold gas into stars, the spectrophotometric evolution of the resulting stellar 
populations, the feedback from star formation and evolution on the properties of prestellar 
gas, and the build-up of large galaxies by mergers. 

As in the case of the CMBR, an arbitrary pattern of fluctuations can be mathematically 
described by an infinite sum of independent waves, each with its characteristic wavelength 
A or comoving wave number k and its amplitude Sk- The sum can be formally expressed 
as a Fourier expansion for the density contrast at comoving spatial coordinate r and world 
time t, 

6{r,t)^J2^k{t)e^-\ (25) 

where k is the wave vector. 

Analogously to Eq. (23) a density fluctuation can be expressed in terms of the dimen- 
sionless mass autocorrelation function 

e(r) = (5(r05(r+rO) oc Ed^^'W H^*'^ • (26) 

which measures the correlation between the density contrasts at two points r and ri. The 
powers deflne the power spectrum of the rms mass fluctuations. 
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Figure 2: Top panel: a compilation of recent CMB data |38]. The solid line shows the result of 
a maximum-likelihood fit to the power spectrum allowing for calibration and beam uncertainty 
errors in addition to intrinsic errors. Bottom panel: the solid line is as above, the solid squares 
1 38 1 and the crosses [^] give the points at which the amplitude of the power spectrum was 



estimated. For details, see reference [pB ]. 
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P{k) = {Mt)\') . (27) 

Thus the autocorrelation function ^(r) is the Fourier transform of the power spectrum. 
This is similar to the situation in the context of CMB anisotropics where the waves 
represented temperature fluctuations on the surface of the surrounding sky, and the powers 
af were coefficients in the Legendre polynomial expansion Eq. (24). 

With the lack of more accurate knowledge of the power spectrum one assumes for 
simplicity that it is specified by a power law 



P{k) oc r= , (28) 

where is the spectral index of scalar fiuctuations. Primordial gravitational fiuctuations 
are expected to have an equal amplitude on all scales. Infiationary models also predict 
that the power spectrum of matter fiuctuations is almost scale-invariant as the fiuctuations 
cross the Hubble radius. This is the Harrison-Zel'dovich spectrum, for which = 1 (n^ 
= would correspond to white noise). 

Since fiuctuations in the matter distribution has the same primordial cause as CMBR 
fiuctuations, we can get some general information from CMBR. There, increasing Ug will 
raise the angular spectrum at large values of i with respect to low i. Support for £ fa 1.0 
come from all the available analyses: combining the results of references , , |Q by 
the averaging prescription in Section 4, we find 



= 1.02 ±0.06 ±0.05 . (29) 

Phenomenological models of density fiuctuations can be specified by the amplitudes 
6k of the autocorrelation function ^(r). In particular, if the fiuctuations are Gaussian, 
they are completely specified by the power spectrum P{k). The models can then be com- 
pared to the real distribution of galaxies and galaxy clusters, and the phenomenological 
parameters determined. 

As we noted in Section 4, there are several joint compilations of CMBR power spectra 



and LSS power spectra of which we are interested in the three largest ones |3q|, |41 



Combining their results for by the averaging prescription in Section 4, we find 



nr. 



0.29 ±0.05 ±0.04 



(30) 



If the Universe is spatially fiat so that i^o = 1? this gives immediately the value Qx = 
0.71 with slightly better precision than above. To check this assumption we can quote 
reference from their Table 5 where they use all data. 



1.00 ± 



0.03 
0.02 



(31) 
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Note, however, that this result has been obtained by marginahzing over all other param- 
eters, thus its small statistical errors are conditional on n^, Qm, being anything, and 
we have no prescription for estimating a systematic error. 

A value for Qx can be found by adding Qx—Qm in Eq- (22) to Qm, thus = 0.79±0.12. 
A better route appears to be to combine Eqs. (30) and (31) to give 

= 0.71 ±0.07 . (32) 

Still a third route is to add Qq and Qx — ^m, or to subtract them, respectively. Then 
one obtains 



0.25 ±0.05 



0.75 ±0.05 



The routes making use of Qx — from Eq. (22) are, however, making multiple use of 
the supernova information, so we discard them. 

Before ending this Section, we can quote values also for wx and go- The notation here 
implies that wx is taken as the equation of state of a quintessence component, so that its 
value could be i^a > — 1- The equation of state of a cosmomological constant component 
is of course wx = —I- In a flat universe wx is completely correlated to f^A and therefore 
also to Qm- 

We choose to quote the analysis by Bean and Melchiorri who combine CMBR 
power spectra from COBE-DMR 0, MAXIMA BOOMERANG DASI §, the 
supernova data from HSST [l^ and SCP [|lTl, the HST Hubble constant quoted in 
Eq. (15), the baryonic density parameter h"^ = 0.020 ±0.005 and some LSS information 
from local cluster abundances. They then obtain likelihood contours in the wx , space 
from which they quote the la bound wx < —0.85. If we permit ourselves to restrict their 
confidence range further by using our value Q^n = 0.29 ± 0.06 from Eq. (30), the result is 
changed only slightly to 



Wx < -0.86 , {la CL) 



(33). 



Finally, the deceleration parameter is not an independent quantity, it can be calculated 
from 



go = ^^m -^x = l^m ~^o = -0.56 ± 0.04 . (34) 

The error is so small because the Qm and the fix errors are completely anticorrelated. 
Note that the negative value implies that the expansion of the Universe is accelerating. 
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Parameters Values References 



Ho 


73 ±7 




i 






0.019 ±0.002 


our 


compilation 




0.29 ±0.06 


our 


compilation 


fix 


0.71 ±0.07 


our 


compilation 


no 


1 n +0-03 

-■-•U -0.02 








Us 


1.02 ±0.08 


our 


compilation 




< - 0.86 


our 


compilation 


Qo 


- 0.56 ±0.04 


our 


compilation 



Table 2: Best values of the dynamical parameters. The errors include la statistical errors and 
our estimates of systematic errors, except for Qq which is statistical only. The Hubble constant 
Ho is given in units of km s^^Mpc^^ 

7 SUMMARY 

Information on the dynamical parameters of the Universe are coming from the Big Bang 
nucleosynthesis, from the fluctuations in the temperature and polarization of the cos- 
mic microwave background radiation, from the large scale structures of galaxies, from 
supernova observations and from many other cosmological effects that may not yet be 
of interesting precision. The results of different analyses are now converging towards 
agreement when in the past disagreements of the order of 100% have been known. 

In this review we have taken the attitude that remaining disagreements reflect sys- 
tematic errors coming either from the observations or from differences in the methods 
of analysis. We have then compiled the most precise parameter values, combined them 
and added our estimates of such systematic errors. This we have done for the baryonic 
density parameter Qb h'^i the density parameter of the matter component VLm, the density 
parameter of the cosmological constant fi^, the spectral index of scalar fluctuations ra^, 
the equation of state of the cosmological constant wa, and the deceleration parameter go- 
In addition we quote the best values of the Hubble parameter Ho and the total density 
parameter VLq from other sources. In Table 2 we summarize our results. 

The conclusion is not new: that the Universe is spatially flat, that some 25% of grav- 
itating matter is dark and unknown, and that some 70% of the total energy content is 
dark, possibly in the form of a cosmological constant. 

Acknowledgements: S. M. H. is indebted to the Magnus Ehrnrooth Foundation for support. 
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